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Abstract 


The purpose of this research is to present a matrix method for solving 
system of linear Fredholm integro-differential equations(FIDEs) of the second 
kind on unbounded domain with degenerate kernels in terms of generalized 
Laguerre polynomials(GLPs). The method is based on the approximation of 
the truncated generalized Laguerre series. Then the system of (FIDEs) along 
with initial conditions are transformed into the matrix equations, which cor- 
responds to a system of linear algebraic equations with the unknown gen- 
eralized Laguerre coefficients. Combining these matrix equations and then 
solving the system yields the generalized Laguerre coefficients of the solution 
function. In addition, several numerical examples are given to demonstrate 
the validity, efficiency and applicability of the technique. 


Keywords: Systems of linear Fredholm integro-differential equations; Un- 
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1 Introduction 


The main object of this paper is to approximate the solution system of Fred- 
holm integro-differential equations of the second kind on a semi-infinite do- 
main of the following form: 


U'(x) = F(x) + an w(t)K(a2,t)U(t)dt, «ER, (1) 


along with initial condition U(0) = A, where p € R, and 
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A= feats Cisterwadiotlingl ° § (2) 


In system (1), w(t) = t*e~*(a@ > —1) and K(za,t) a function of two vari- 
ables x and t, is called the kernel that might have singularity in the region 
D = {(a,t):0 < x,t < co} and F(z) is continuous function and A is fixed 
constant vector, and U(x) is the unknown vector function of the solution that 
will be determined. The considered equation arises in a number of important 
problem of elasticity theory, neutron transport, particle scattering and the 
theory of mixed-type equations [11,13,17]. System of linear Fredholm integro- 
differential equations of the second kind on unbounded domain can not be 
analytically solved easily. Therefore, it is required to obtain the approxi- 
mate solutions. It’s the reason of great interest for solving these equations. 
But numerical methods includes Quadrature, Petrov-Galerkin, Nystrom and 
Galerkin methods with Laguerre polynomial as a bases function for solving 
infinite boundary integral and integro-differential equations are used ago that 
their analysis may be found in [1,7,9, 10, 12,16]. On the other hand, there 
are several numerical techniques for solving fractional differential equations 
(FDEs) on the half line using generalized Laguerre polynomials [2-6]. How- 
ever, method of solution for equation (1) is too rear in the literature. In 
the present work, we are going to use the operational matrix of generalized 
Laguerre polynomials to find the approximate solutions for the system of 
FIDEs on the half-line. Next sections of this paper are organized as follows: 
In Section 2, we describe some necessary definitions and give some relevant 
properties of the GPLs which is required for our subsequent development. 
Section 3, is devoted to the approximation of the function f(x) and also the 
kernel function k(x,t) by using GPLs basis. Also the upper bound of the 
approximation error is presented. In Section 4, we obtain the operational 
matrix of integration by GPLs. In Section 5, we implemented the matrix 
method on the system of linear Fredholm integral-differential equations on 
unbounded domain and convert them to a linear algebraic system of equa- 
tions. In Section 6, presented numerical examples that shows the efficiency 
and accuracy of the proposed method. Also a tall conclusion is given in 
Section 7. 


2 Preliminaries and basic definitions 


In this part, for the reader’s convenience, we give some basic definitions and 
properties of the generalized Laguerre polynomials, which are used further in 
this article. 
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Let Ry := A = (0,00) and w(x) = x%e-* be a weight function on A in the 
usual sense. We define the following: 
L? 


“(a (A) ={v:v is measurable on A and |lv||y a < co}, (3) 


equipped with the following inner product and norm: 
dl 
(U,V) wa = i u(a)o(x)w (x)dx,  |lvlluco = (¥,¥)2ca)- (4) 
A 


Next, suppose L& (a) be the generalized Laguerre polynomials of degree n, 
defined by the following: 


1 


L‘*) (x) = nite Or (et) n= Oia (5) 


L‘) (x) (generalized Laguerre polynomials) are the nth eigenfunction of the 
Sturm-Liouville problem: 


! / 
aXe” (ene (L2@)) ) + AnL (x2) =0, «eA, (6) 


with the eigenvalues A, = n [8, 14]. 
Generalized Laguerre polynomials are orthogonal in L?,.,)(A) Hilbert space 


an,-x 


with the weight function w'~) (x) = a%e~* satisfy in the following relation 
[ xe" Ly” (x) LO) (x)dz = Yn n,m; Vn,m = 0, (7) 
0 


a — T(n+a+1) 


n= Torti) The explicit 


where dy, is the Kronecher delta function and 7, 
form of these polynomials is in the form 


LY (a) = y Eee", (8) 
i=0 
where 
n+a : 
. | (—1)' 
pent las (9) 


i! 


These polynomials are satisfied in the following recurrence formula 


LM (e)=1, LO (2) =1+a-z, 


[(2n ta+1—2)L'%) (2) — (n+ a)L&©, (x) sb Le Deer: 
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The case a = 0 leads to the classical Laguerre polynomials, which are used 
most frequently in practice and will simply be denoted by L,,(a). An im- 
portant property of the Laguerre polynomials is the following derivative re- 
lation [10]: 


n-1 


(29a) = 1a). (11) 


i=0 
(a) (¥) F F 
Further, ( L; (x)) are orthogonal with respect to the weight function 
weet) (x), ie. 


| (LP) (2) (LP) (2) (a)de = 92465; 5, Vij >0, (12) 
0 


where 7°*? is defined in (7). 


n 


3 Approximation of functions by using GLPs 


An arbitrary function f(x) € L?(.)(A) may be expanded into generalized 
Laguerre polynomials as: 


Co 


f(z) => OL (a), (13) 


i=0 


where the generalized Laguerre coefficients FO are given by 


co (a) a,x 
oo L;"’ (2) ae — 
ie =| (48) Tat) f(a)dx, ti=0,1,.... (14) 


u 


The series converges in the associated Hilbert space L? ... (A), iff 


= CO pea ~ ©9 ite (a) 
Wilts PO ep e@tdr= LLM cee.) 


In practice, only the first (n + 1) terms of generalized Laguerre polynomials 
are considered. Then we have 


f(a) & SS AOL! (@) = FTL, (16) 
i=0 


where the generalized Laguerre coefficient vector F' and generalized Laguerre 
vector L, are given by as follows: 
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Fa [fA s es OVP, Le = L57(@), L(@)s---, LOI? 7) 


Now in the following lemma we present an upper bound to estimate the error. 


Lemma 1. Suppose that the function f : A > R isn+1 times continuously 
differentiable (i.e. f € C"**(A)), andY = Span{L\ (2), L(x), miley L (2)}. 
If F™ Lz be the best approximation f out of Y then mean error bound is pre- 
sented as follows: 


Nv/(2n+ a+ 2)! 
(n+ 1)! j 


If -FLellez wy S (18) 


where N = maxzcea pr a) 


Proof. We know that the power basis {1, «,...,2”} forms a basis for the space 
of all polynomials of degree less than or equal to n. Therefore, we define 


yi(x) = f(0) + af’(0) 4 = f"(0) +---4 a f((0). From Taylor expansion 
we have i 
fa) — ni @)l SO (ne) Tah (19) 


where 7, € (0,00). Since F? L, is the best approximation f out of Y, y, € Y 
and using (19) we have 


N?(2n +a + 2)! 


Tyr 2 2 
|f-F Lalln2 (a) Se = villae cay < GE=E (20) 
Then, by taking square roots we have the above bound. 
This Lemma shows that the error vanishes as n — oo. 
We can also approximate the function of two variables, k(,t) € L? ..)(A’) 
as follows: — 
k(a,t) & S> S71 (wy MOL (t) = LEK Ly. (21) 


, for i,f=0,1,...,n, (22) 


so that, (.,.) denotes the inner product. 
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4 Operational matrix of integration, development and 
applications 


The main objective of this part is to obtain the operational matrix of the 
integration by GPLs. 


Theorem 1. Suppose L, be the generalized Laguerre vector defined in (17) 
then, 


/ L,dt ~ PL, (23) 
0 


where P is the (n+1) x (n+1) operational matrix for integration as follows: 


Q(0,0,a) (0, 1,a) Q(0,2,a) +--+ O(0,n, a) 
Q(1,0,a) Q(1,1,a) Q0,2,a) --- AC, n, a) 
Pl Oe One Lae ais Ema | 4) 
OGn, 08) CG.) Opa = Oe) 
where 
ne es ae (-1)**"jTG+ta+)r(k+at+r42) 
Gj¢)= >) 2 GomiG =n DT EKO Te coda 


Proof. (See [4]). 


Also, we can see the extent of this theorem for solving fractional differen- 
tial equations. For example, see [2, 5, 6]. 


5 Implementation of the matrix method 


In this section, we solve the system of linear Fredholm integro-differential 
equations of the second kind on unbounded domain (1). To this end, we 
consider the ith equation of (1) as follows: 


u(x) = file)+e | te! S > hig (a, t) uj (t)dt, us(0) =a;,i=1,...,m, (26) 
j=l 


where f; € L?(.)(A), kij € L?(.)(A”), and u/(x) represents the first order 
. . wy . ey . . eye 
derivative of u;(x) with respect to x, a; are constants that give the initial con- 
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ditions and u,; is an unknown function. In order to approximate the solution 
of equation (26), we approximate functions f(a), u;(v) and k,;(x,t) with 
respect to generalized Laguerre polynomials as mentioned in the previous 
section as follows: 


filz) 2 FPL,, ui(c) 2 CP Ly, u(0) 2% CRLe, kiy(z,t) ~ LE KijL,, (27) 


where F;, Ci for i = 1,...,m are known (n+ 1) x 1 vectors and Kj; for 
i,j =1,2,...,m are known (n+ 1) x (n+ 1) matrices. Then, for i = 1,...,m, 
we have: 


u(x) = | ul (t)dt + u;(0) ~ | OF Trdt + C2 Ly» (CP P+C4) Le; (28) 


where P is a (n+ 1) x (n+1) operational matrix of integration given in (23). 
By substituting the approximations (27) and (28) into equation (26), we get 
the following: 


Cai Rk + a tet S 0 LE KizLtLt (p" Ch + Cjo)dt 
0) 


j=l 


= LPF, + pL? S_ Kis | te" Ly Lf dt} (p™ Cy + Cjo) 
j=l 


= LT F, + pL? S~ KijQ(p" Ci + Cjo). (29) 


j=l 


Then we have following system of linear equations: 


ClL=F.+ p>) KiyjQ(P™C)+Cjo), 4=1,...,m, (30) 
j=l 
where ss 
Q= | te LL dt = (a), i,j =0,1,...,n, (31) 
0 
and Q is a (n+ 1) x (n+ 1) matrix with elements 
gi = 7 eT (tL (at, 4, f = 0,1,..., 0. (32) 
0 


By solving the linear system of algebraic equations (30), we can achieve the 
vector Ci for i =1,...,m, then we will have 


CF HCP PSC) SU aeci. t= Tekin (33) 


That are the approximate solution for our system of (1). 
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6 Numerical Examples 


In this section, we give several illustrative examples for demonstrate the 
efficiency of our proposed method to approximate the solutions system of 
Fredholm integro-differential equations of the second kind along with initial 
condition on a semi-infinite domain. For each example, we find the approx- 
imate solutions using different degree of generalized Laguerre polynomials. 
The results obtained by the present method reveal that the proposed method 
is very effective and convenient for system (1) on the half line. In all exam- 
ples, the package of Matlab (2013) has been used to solve the test problems 
considered in this paper. 


Example 1. For the first example, consider the following of Fredholm 
integral-differential equation on unbounded domain (constructed): 


247131410303000045 »  38903199231847830919 


/ —_ 
(2) = ~36998797018063068 ° ~~ 144115188075855872 


+| tze-*(a? + t?)u(t)dt, u(0) =1. (34) 


Exact solution of this problem is u(a) = 23 —22+1. If we apply the tech- 
nique described in the section 5, with a = 4 and n = 3, then the approximate 
solution can be expanded as follows: 


3 


u(x) ~ Sof? {9 (2) = CT Le, (35) 
i=0 
where 
C= [h™, of, cf, SO, (36) 


Hence, from Eqs. (16), (21), (23), and (31), we find the matrices 


—125363/424 15/2 —5 20 3/2 -10 0 
| 17011/496 _| -5 000 _| 3/8 1-10 
ne —6434/469 cea Wy ae (cae 5/16 0 1-1]? 

0 0 000 35/128 0 0 1 

148/167 0 0 0 ] 

Q= 0 222/167 0 0 
~ 0 0 555/334 0 
0 0 0 2053/1059 


Next, we substitute these matrices into equation (30) and then simplify to 
obtain 
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one —119/5475 —19/7262 —574/2251 —865/5409] [ —13873/48 
cl) | _ | 161/2349 347/1578 93/632 —247/3501| | 4211/141 
(oy | ~ | —181/6602 552/1769 1487/1580 193/6839 |" | —6427/538 
2 
(a)! 0 0 0 1 0 
C3 

(37) 


By solving the linear system (37), we have the following: 
a Si 37 a v a t a ti 
oy =, a = 15, ey” =6, of =0. (38) 


By substituting the obtained coefficients in (33) the solution of (34) becomes 


89 (a 97 (a i a 
u(r) ~ Lh” (x) — LY" (a) + 21Ly" (2) — 615” (2), (39) 
or briefly 
u(x) ~ 2 — Qe +1, (40) 


which is the exact solution. Also, if we choose n > 4, we get the same ap- 
proximate solution as obtained in equation (40). Numerical results will not 
be presented since the exact solution is obtained. 


Example 2. As the second example, we consider the following system of 
linear Fredholm integro-differential equations on unbounded domain (con- 
structed): 


ui (x) = fi(x) + i t2e~* (Qa + #7)(uz(t) + u2(t))dt, 


Co 
1 
u(x) = fo(x) +f t2e~*(t — a*)(ur(t) — ua(t))dt, (41) 
0 
_ 9.2 _ 87307746120759955 6631788499575074881 
where f(x) = 3x DBI TOISISGSS248 © — “TBolas9sn0948T984 and 


f (x) = 98782478468059837 ,.2 L Dp 853121404951425865 
2 9007199254740992 ' 18014398509481984 * 


Subject to the initial conditions u (0) = 1 and u2(0) = 1. The exact 
solutions of this problem are ui(x) = x? + 2x +1 and u2(x) = 2? +1. If 
we apply the technique described in this paper and solve equation (41) with 
a = 4 and n=3. For this system we get: 


137 7 (a 113 (a a a 

uy (x) = bs (x) — wae. M(x) + 21L6) (x) — 6L&) (x) = 23 + 22 +1, 
19 a a a a 

g(a) = Ly (e) — SLY (@) + 204 (#2) + (LS(e) =a? +1, (42) 


which is the exact solution. Also, if we choose n > 4, we get the same ap- 
proximate solution as obtained in equation (42). Numerical results will not 
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be presented since the exact solution is obtained. 


Example 3. As the third example, consider the following system of linear 
Fredholm integro-differential equations on unbounded domain (constructed): 


ui (a) = fila) + ie te *~*(sin(t — x)u;(t) + tug(t))dt, 


U(x) = fo(z) +f te~" (xtu,(t) — e~*ug(t))dt, (43) 
0 
with fi(~) = 1— $(1+ 2sina + 2cosx)e~* and f2(x) = —6x — $e~* and 
with the exact solutions uj(”) = x, u2(%) = e~* and boundary conditions 
u1(0) = 0 and u2(0) = 1. We apply the generalized Laguerre series ap- 
proach and solve equation (43). Table 1 shows the absolute values of error 
\e| = |u2 (a) — Ua (a)|, where u(x) is the exact solution of equation (43) and 
Ug(x) is the approximate of uo(x#) for n = 20, and n = 30 with a = 1 using 
the described method in equally divided interval [0,1]. Note that absolute 


Table 1: Absolute errors for Example 3 


a Ly n = 20 n = 30 

0 | 0.0 | 5.4836e — 006 | 7.6834e — 009 
1 | 0.1 | 1.2376e — 006 | 3.6294e — 010 
2 | 0.2 | 4.7027e — 007 | 1.1165e — 009 
3 | 0.3 | 8.2669e — 007 | 5.9280e — 010 
4 | 0.4 | 5.7470e — 007 | 1.8341e — 010 
5 | 0.5 | 1.491le — 007 | 5.9768e — 010 
6 | 0.6 | 2.2241e — 007 | 5.8849e — 010 
7 | 0.7 | 4.4523e — 007 | 3.0976e — 010 
8 | 0.8 | 5.0477e — 007 | 4.8167e — 011 
9 | 0.9 | 4.2952e — 007 | 3.3937e — 010 
10 | 1.0 | 2.6633e — 007 | 4.8646e — 010 


errors for ui(a) is zero. 


Corollary: If the exact solution of the system (1) be a polynomial, then the 


proposed method will obtain the real solution. 


Example 4. 


u(r) =e * — ive+ | t2e~*V/axtu(t)dt, u(0) =1. 


oo 
0 


Our last example is following of linear Fredholm integro- 
differential equation on a semi infinite interval (constructed): 
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With the exact solution u(x) = 2—e7*. In Table 2, the numerical results 
of the presented method at some selected nods for n = 10, and n = 12 are 


displayed. 


Table 2: Absolute errors for Example 4 


a Ls n= 10 n=12 

0 | 0.0 | 1.3000e — 003 | 3.6116e — 004 
1 | 0.1 | 4.4620e — 004 | 9.3426e — 005 
2 | 0.2 | 7.1845e — 005 | 4.5500e — 005 
3 | 0.3 | 3.2705e — 004 | 9.9666e — 005 
4 | 0.4 | 4.0407e — 004 | 1.0182e — 004 
5 | 0.5 | 3.6812e — 004 | 7.5666e — 005 
6 | 0.6 | 2.6834e — 004 | 3.7765e — 005 
7 | 0.7 | 1.4076e — 004 | 8.8208e — 007 
8 | 0.8 | 1.0787e — 005 | 3.3527e — 005 
9 | 0.9 | 1.0468e — 004 | 5.6602e — 005 
10 | 1.0 | 1.9540e — 004 | 6.8843e — 005 


7 Conclusion 


Obtaining the analytic solutions for system of linear Fredholm integro- 
differential equations of the second kind, along with initial conditions on 
unbounded domain are usually difficult. In many cases, it is required to ap- 
proximate solutions. For this reason, a new matrix approach which is based 
on the generalized Laguerre operational matrix of integration is proposed. 
The solution procedure is very simple by means of generalized Laguerre poly- 
nomials expansion and only in a few terms lead to high accurate solutions. 
The main goal of the presented technique was deriving an approximation to 
the solution system of linear Fredholm integro-differential equations on un- 
bounded domain. To illustrate the method and its efficiency, four examples 
were provided. In the first and second examples, we obtained the exact so- 
lution. Another considerable advantage of the method is that the nth-order 
approximation gives the exact solution when the solution is polynomial of 
degree equal to or less than n. If the solution is not polynomial, generalized 
Laguerre series approximation converges to the exact solution as n increases. 
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